Quantized squeezing and even-odd asymmetry of trapped bosons 
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We investigate the exact nature of the superfluid-to-Mott-insulator crossover for interacting 
bosons on an optical lattice in a one-dimensional, harmonic trap by high-precision density-matrix 
renormalizat ion-group calculations. The results reveal an intermediate regime characterized by a 
cascade of microscopic steps. These arise as a consequence of individual boson "squeezing" events 
and display an even-odd alternation dependent on the trap symmetry. We discuss the experimental 
observation of this behavior, which is generic in an external trapping potential. 
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The recent rapid developments in ultra-cold-atom ex- 
periments for both bosons and fermions have greatly 
stimulated the fuller exploration of a number of fun- 
damental properties of strongly correlated systems p. 
One of the most important phenomena to be observed 
and characterized is the transition from a superfluid to 
a Mott-insulator phase displayed by interacting bosonic 
atoms in optical lattices [2]. The Bose-Hubbard model 
[3] shares a number of the properties of the Hubbard 
model for correlated electrons, whose properties, in- 
cluding superconductivity and Mott-insulating behavior, 
have challenged condensed matter physicists for half a 
century. This model can be realized in one, two, and 
three dimensions on an optical lattice [4 , with the enor- 
mous advantage over electronic solids that the ratio be- 
tween kinetic and interaction parameters of the particles 
is in principle continuously tunable. 

In experiments, a trapping potential is required to con- 
fine the atoms. It was shown both theoretically |5 and 
experimentally [11 El |7] that the pure Mott-insulating 
phase, where the boson distribution is uniform despite 
the trap profile, is obtained only very at high interac- 
tion strengths. Otherwise, a spatial "shell structure" of 
the local density is found [11 [HI [9] , with Mott-insulating 
and superfluid regions present simultaneously in regions 
of different trap depth. Extensive studies of this behav- 
ior include a demonstrated compression of the superfluid 
region as the Mott insulator is approached JO] and quan- 
titative calculations both of the shell structure [Til [E] 
and of accompanying features in the visibility [T3] . 

Here we show that the nature of the Mott transition in 
a trapping potential is considerably more complex than 
recognized previously: "superfluid" bosons are expelled 
one by one into the Mott-insulating regions with in- 
creasing interaction strength. The intermediate regime 
where this sequence of microscopic quantum transitions 
takes place is the necessary consequence of the additional 
energy scale introduced by the trap, and the resulting 
physics is qualitatively different from the infinite system. 

To explore the nature of the states intermediate be- 
tween the perfect superfluid and the Mott insulator in 



a trap, we consider the one-dimensional (ID) Bose- 
Hubbard model with a harmonic trapping potential. This 
system is realized experimentally for a cloud of interact- 
ing bosons loaded in an optical lattice with one dominant 
recoil energy. The Hamiltonian is expressed [3l |4] as 

H = -ty^ f^l^i+i + ^?7y^ni(ni - 1) 

i=l,L i=l,L 

+VTj2ni[i-iL + l)/2f, (1) 

i=l,L 

where the lattice spacing is set to 1 and the on-site in- 
teraction U and (bond-centered) trapping potential Vt 
are measured in units of the hopping integral t. We ap- 
ply the density-matrix renormalization-group (DMRG) 
technique [14 , using a specially modified finite-system 
algorithm optimized for the accuracy of the low-lying 
energy states. An appropriate sweeping procedure is es- 
sential to guarantee convergence in this delicate problem 
^15j . Sweeping is conducted gradually from the middle 
to the two ends, convergence is obtained for each sweep- 
ing length / G [4,L], and the three lowest eigenstates are 
targeted simultaneously. Extremely high accuracies are 
required to allocate reliably the symmetries of degenerate 
states, and these were examined [16] for different system 
sizes L, numbers n of bare states per site, and numbers 
m of states retained per block, for each total boson num- 
ber N and interaction U. Working up to n = 32 and 
m = 1200, we found that n = 8 and m = 100 are suf- 
ficiently large for G [4, 14] to ensure truncation errors 
smaller than the symbol sizes in the figures below, while 
n = 32 and m = 400 are necessary for /7 < 4. 

Systematic studies were performed for systems with 
20 < TV < 70 bosons in a trap with Vt = 0.01, for which 
L = 80 is sufficient to remove any boundary effects. The 
unconventional properties of the intermediate regime be- 
gin at A/" = 40 and become progressively clearer with 
increasing N until a higher Mott plateau is reached. We 
have performed many calculations to amplify the current 
results [H], but this parameter and filling regime, shown 
in Fig. 1, ensures a minimal system reflecting all of the 
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FIG. 1: (color online). Results for Vt — 0.01, t = 1, and system length L = 80 for AT = 30, 40, 50, and 60 particles, (a) Density 
profiles rii as functions of U . (b) Energy gaps Ai,2; inset: Ai,2 for low U . (c) Visibility v (left axis) and peak width w (right). 



qualitative physics. We focus only on the ground state 
and not on the nature of the excitations [16 ; however, 
because /7-induced changes in the ground state appear as 
(avoided) level-crossings, it is instructive to consider the 
lowest three energy levels, £^0,1,2, represented in Fig. 1 by 
the gaps A 1^2 = ^1,2 — Eq- We stress again that this is a 
finite system, both theoretically and in experiment: we 
provide numerically exact results for the intrinsic physics 
of this system, and there are no "finite-size effects" in 
the conventional sense of approximations to a thermo- 
dynamic limit. We also calculate the particle numbers 
(n^'^'*) for single, double, and triple occupation, the vis- 
ibility z^, and the peak height ^max and width w (taken 
as the full width at half the sum of the maximum and 
minimum values) of the momentum distribution function. 
The measurable physical quantities v and w [1 illustrate 
most clearly the effects we consider here. 

Figure 1 shows three clear regimes. At low U is the 
bell-shaped distribution of a "superfluid" (SF) system, 
which has strong visibility and a trap-related gap above 
the ground-state condensate. At high U is the fiat dis- 
tribution {rii = 1), despite the energy cost of the trap at 
the edges, of a pure "Mott insulator" (MI) regime, with 
falling visibility and gaps of order t to states at the edges. 
Between these two limits is a complex intermediate phase 
characterized by oscillating and vanishing energy gaps, 
step-like changes in the particle distribution and the vis- 
ibility, and a spatial "shell structure" with MI plateaus 
around a central SF region. Note that we use the abbre- 
viations SF and MI to refer not to the bulk phases but 
respectively to the local regions of the distribution shown 
by the red-yellow curve and by the green plateaus. The 
almost constant shape of the blue "tails" in Ui at the 
system edges is dictated not by U but by Vt- 



For the homogenous Bose-Hubbard model, the Mott 
transition occurs at Uc = 3.61±0.13 [5 . In a trap, there 
is little evidence of the energy scale Uc, although it cor- 
responds loosely to the onset of oscillations in Ai^2 and 
in Hi. It is easy to deduce that a pure MI phase is es- 
tablished only when the energy cost U of the last boson 
in the central SF region exceeds that of placing it at 
the edge. One effective description of this process is to 
write Uc2 - a{l)t 0^ IVt{N - 5)^, where a{l) - 4 [T6| 
represents a kinetic-energy contribution due to the last 
boson and 6^6 [Fig. 2(b)] represents the edge effects at 
lowest order. Both a and 6 are corrections to a leading 
quadratic dependence of Uc2 on A^, arising directly from 
the harmonic trapping potential, which makes the width 
in U of the intermediate regime very significant. This 
illustrates directly how Vt acts as an additional energy 
scale, which is responsible for the presence of the addi- 
tional, intermediate regime. Clearly, if Vt = the atoms 
are no longer trapped, or their density controlled: the 
trap (Vt) is a non-pert urbative term in ([l]), ensuring 
that any ensemble of trapped atoms is generically in a 
different class from the homogeneous, infinite system. 

As shown in Fig. 1, the width in U of the intermediate 
regime indeed increases strongly with N. The oscilla- 
tions in A 1^2 and rii grow with U until level-crossings 
and steps appear (we refer to "plateaus" in i and "steps" 
in U). In this cascade of steps, Ai vanishes over short 
but finite ranges of [/, indicating the double degeneracy 
of the ground state when the two MI plateaus contain an 
odd number of particles. Each level-crossing represents 
a squeezing process where one particle is ejected from the 
central SF region and placed on the outer edge of one or 
other plateau: these processes are quantized. The loca- 
tions and number of these crossings, which correspond 
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FIG. 3: (color online) Lowest energy gaps for = 60, L = 80, 
and Vt — 0.01 with trap center offsets 0.25 (a) and 0.5 (b). 
The case with offset is shown in Fig. 1 (center right panel). 



FIG. 2: (color online). Superposed distributions of particles 
in (a) doubly occupied sites (nf ) and (b) singly occupied sites 
(n- ) for iV = 60, L = 80, and Vt = 0.01. Shown in each panel 
are 47 curves for values 7.6 < U < 12.2, which fall into 10 
different sets corresponding to 0, ... 9 bosons in the central SF 
region. Red curves for two remaining SF bosons and green for 
one. (c) Schematic representation of generic "state diagram" 
for a trapped boson system as a function of U and N. Solid 
hues mark finite-system "phase transitions," dot-dashed hues 
a typical filling Na with two n=l Mott plateaus. 

directly to jumps in , and w as functions of [/, are 
not universal and nor are their spacings identical: the 
examples in Fig. 1 favor even numbers of particles in the 
central SF region ("even-odd effect"). We note that such 
series of discrete, single-particle transitions are present in 
Monte Carlo simulations performed for a harmonic trap 
with n = 1 and n = 2 MI plateaus [13] and in DMRG 
calculations for a double-well trap [17]; however, neither 
set of authors appeared to recognize their significance as 
quantized squeezing, or noticed an even-odd effect. 

Additional insight into the squeezing process can be 
found in Figs. 2(a) and (b), which show the evolution 
with U of Ui across the intermediate regime for = 60 
(right panels of Fig. 1). While nj and higher site oc- 
cupations are always negligible, n| and nf are perfectly 
anticorrelated. A key qualitative point must be made 
here: although = 1 in the MI plateaus, this is not 
due to perfect single-site occupation (because nf 7^ 1), 
and the finite constant nf reflects the complete coher- 
ence of the many-body wave function mediated by con- 
tinuous kinetic processes coupling the central SF to the 
edges through the MI regions. In spite of this system- 
wide coherence, does not evolve continuously, but in- 
stead undergoes 9 abrupt changes across the intermediate 
regime [Figs. 2(a) and (b)] as U acts to squeeze the ini- 
tial 9 bosons in the central SF region successively down 
to zero. The areas represent the exact boson numbers, 
giving complete quantitative information on the spatial 
extent and the structure of the central SF part of the 
wave function, the number of weak maxima for 2, 3, . . . 
particles reflecting the boson phase coherence. 

Returning to the quantitative information in Fig. 1(b), 
and to the question of a "phase diagram" for the fi- 
nite system of trapped bosons, we define the onset Ud 



of an intermediate phase, which corresponds to the ap- 
pearance of two true MI plateaus as opposed to flatten- 
ing shoulders, from the first level-crossing. Proceeding 
as above, U,, - a(iVl, + l)t IVriN - 5 - N^,J^ 
where A^^ax number of bosons accommodated in 

the central SF region at Ud and —a{N^^^-\-l)t represents 
the kinetic energy of the squeezed boson. The leading 
dependence of the width of the intermediate regime is 
then Uc2 — Uci ~ ^^rA'A'max^ shown schematically 
in Fig. 2(c). This type of "canonical phase diagram," 
or "state diagram" [11 , is appropriate for cold-atom ex- 
periments where a fixed number of particles is loaded 
and the optical lattice tuned adiabatically. Figure 2(c) 
also shows that no MI can be formed if A" is too low 
(A' < A"^), while for A^ sufficiently large, higher plateaus 
will appear. Although the functional forms and relative 
positions of the phase boundaries can be changed by al- 
tering the trap shape, this is the generic state diagram 
of any ensemble of bosons whose spatial distribution is 
contained by a finite external potential. 

A further essential qualitative feature of single-boson 
squeezing is its even-odd asymmetry. Figure 1 shows 
even boson numbers A" in a bond-centered ("even") trap, 
where the steps of even boson occupation in the central 
SF region are wider (more stable) than the odd ones. 
However, if A" is odd, the situation is exactly reversed, 
and the odd-occuption SF steps are wider [16] . To inves- 
tigate this effect systematically, we move the trap mid- 
point continuously from the bond-centered position to 
a site-centered one ("odd trap"). Again the situation 
is reversed in the site-centered limit (offset 0.5), with 
odd steps more stable [Fig. 3(b)] for even A". This sug- 
gests a very straightforward interpretation: as U is in- 
creased, two states, with the squeezed boson at one side 
of the system or the other, simultaneously become more 
favorable than the states at the center. This picture is 
exact in the t ^ limit, where it is clear analytically 
that bosons are squeezed out of the central SF region in 
pairs, while the effect of a finite kinetic-energy term t 
is to stabilize a narrow (in U) step of the intermediate, 
symmetry-disfavored boson number. Whether even or 
odd steps are favored is a simple consequence of the trap 
symmetry and filling. The confirmation that for spinless 
bosons there are no special, pairwise correlation effects 
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FIG. 4: (color online) Results shown in the format of Fig. 1 
for N — 50 and 60 particles in a 100-site trap as a function of 
the optical lattice depth parameter Vb for Er = 1 (see text). 



in the SF or MI regions is obtained by following the evo- 
lution of the asymmetry with the trap offset: specifically, 
there is no even-odd effect at all for offset 0.25 [Fig. 3(a)], 
confirming its origin in the trapping potential. 

Having understood the physics of quantized squeezing, 
it is possible to comment on its manifestations in a less 
minimal context. The same squeezing effect clearly op- 
erates with more bosons in the trap, within and between 
more complex shell structures of higher MI plateaus. We 
also find an even-odd effect when squeezing bosons from 
central SF or MI regions to the n = 2 or higher Mott 
plateaus, and not only to the system edge (the n = 1 
plateau). Quantized squeezing processes occur in all 
dimensions as individual particles within the coherent, 
many-body wave function are pushed successively out of 
the central SF region with increasing U. In dimensions 
greater than one, much more elaborate even-odd effects 
would be expected: as one example, the generic state de- 
generacy for a square lattice in a circular trap (2D) is 
8, and thus cascades of seven narrow steps might be ex- 
pected between wide 8th steps for a centered trap. For a 
cubic lattice and spherical trap (3D), this number is 48. 
A more challenging question in the physics of strongly 
correlated systems is posed by the situation with finite- 
spin bosons or several boson species in the trap, where 
interaction effects may lead to intrinsic pairing (or more 
complex n-merization) in the SF, or indeed the MI, re- 
gion. This would generate a completely separate mech- 
anism for the tendency of particles to move from SF to 
MI regimes in groups rather than singly. 

We close by discussing how quantized squeezing and 
even-odd asymmetries would be observed in experiment, 
restricting our considerations to ID. Experimental mea- 



surements on optical lattices are performed by varying 
the lattice depth Vb in units of the recoil energy Ej^ . 
Following the calculations of Ref. [7 , we have taken 
t/Er = 1.43(Vo/^r)°-^^exp[-2.07v1yE;], U/Er = 
0.0386(Vb/£;r)°•^^ and Vr/Er = (5.332 + 3.427Vb/£;r) x 
10~^. Clearly, varying Vo/Er explores one curve in the 
2D space of the parameters U and Vr for each TV. Fig- 
ure 4 shows the numerical results obtained as a func- 
tion of Vb for N = 50 and 60. It remains possible, but 
not simple, to find parameters and trap fillings for which 
all three regimes can be observed, although this requires 
scanning over a rather broad range of Vq. While the 
boson-squeezing plateaus are widely spaced in Vb, their 
even-odd asymmetry becomes very pronounced. As a 
final caveat, some concern has been raised recently [18] 
over the heating effects occurring on changing the lattice 
depth for pre-trapped bosons due to energy absorption 
from the optical field, and appropriate measures are re- 
quired to remain sufficiently close to the adiabatic limit. 

In summary, we show that there exists a regime 
intermediate between the pure superfluid and the 
Mott-insulator phase for interacting bosons in a one- 
dimensional optical lattice with a harmonic trapping po- 
tential. We have investigated the sequence of microscopic 
quantum transitions across this regime. These processes 
reflect the quantized squeezing of bosons into the Mott- 
insulator region, and occur as discrete events despite the 
overall coherence of the many-body wave function. The 
even-odd asymmetry of the squeezing phenomenon is a 
direct consequence of the trap symmetry. The "phase 
diagram" of the trapped system is intrinsically, qualita- 
tively different from the infinite system, the intermediate 
regime being the necessary consequence of the additional 
energy scale introduced by the trap. The effects revealed 
by detailed studies of this sort mandate a broader inter- 
pretation of the "Mott transition" in cold-atom systems. 
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